The Minimal Landau Background Gauge on the Lattice 
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We present the first numerical implementation of the minimal Landau background gauge for 
Yang-Mills theory on the lattice. Our approach is a simple generalization of the usual minimal 
Landau gauge and is formulated for general SU(N) gauge group. We also report on preliminary 
tests of the method in the four-dimensional SU(2) case, using different background fields. Our tests 
show that the convergence of the numerical minimization process is comparable to the case of a null 
background. The uniqueness of the minimizing functional employed is briefly discussed. 



In Ref. [1] Cornwall pleaded with the lattice commu- 
nity for an answer to the following question: Can you 
find a way of doing lattice simulations in the background- 
field Feynman gauge ? The reason for this request is that 
one can show [2] — to all orders in perturbation the- 
ory — that there is a simple correspondence between 
the background-field method in the Feynman gauge [3] 
and the so-called pinch technique [4], which allows one 
to build gauge-invariant off-shell Green functions in the 
continuum. 

Let us note that the numerical implementation of Lan- 
dau gauge fixing (e.g. for the evaluation of Green func- 
tions) is well understood [5] . Recently, it has been shown 
that practical simulations of the linear covariant gauge 
are also possible [6] and that, with a suitable discretiza- 
tion of the gluon field, it becomes feasible to treat the 
Feynman gauge [7]. Here we present the first numeri- 
cal implementation of the minimal Landau background 
gauge on the lattice. Our proposal is based on Ref. [8], 
which considers this gauge in the continuum. 

The covariant background gauge condition is intro- 
duced [9] by splitting the (continuum) Yang-Mills field 
A^x) into a quantum fluctuation component Q^x) and 
a background field B^(x), i.e. 



A^x) = Q^x) + B^x) , 



(1) 



where A^ (x) is given in terms of the generators T\> of the 
SU(N) gauge group by A /1 (x) = A b {x)T}, [and similarly 
for Qfj,(x) and B^(x)\. Note that B^(x) is in principle 
arbitrary [10]. Then, the usual covariant gauge condition 



dMx) = \(x) = A b (x)T b 



(2) 



becomes 



dpQ^x) + i[B^x), Q^x)] ee Dp[B\Qp{x) = A(x) 



Here D^[B] is the background- field covariant deriva- 
tive and A b (x) is a Gaussian-distributed real variable. 
Clearly, for a null background field B^(x) = one has 
Qfj,{x) = Afj,[x) and the usual covariant gauge condition 
(2) is recovered. For A(x) — the gauge condition (3) is 
the Landau background gauge condition. 

Let us recall that the continuum gauge transformation 
of the Yang-Mills field, i.e. 

AW(x) = g{x)A, l {x)g\x) - ig(x)d^{x) , (4) 

becomes 



A<*)(x) « A^x) + D^AMx) 
if an infinitesimal gauge transformation 



g(x) = aq>[-ij(x)'_ 



i 7(2;) 



(5) 



(6) 



is considered, where 7(2;) = r ) h (x)T\ } . [Note that, with 
our notation, the generators Tf, are Hermitian. In what 
follows we will also employ the relations TrTJ, = and 
Tr {TfcTc} cx Sbc-] Then, using the splitting in Eq. (1), 
there is clearly no unique way of defining the infinitesi- 



mal gauge transformations (x) and B^'(x) for the 
quantum fluctuation and the background fields. In- 
deed, depending on which of the three terms d fJl j(x), 
i [Q^ix), 7(2;)] and i [B^(x), j(x)] [see Eq. (5)] are in- 
cluded in Q\?\x) and Bjf\x), eight different sets of 
gauge transformations arise naturally. Among these, two 
common choices are 

Q^ ) {x)=Q ll {x) + D^BMx) + 1 [Q„(aO, 7 (aO] (7) 

B^(x)=B u (x) (8) 



tin), 



and 



)(s)(V\ = 



(3) 



(x) = Qn(x) + i [Qn(x), 7(2); 
B^{x)=B^x) +D ll [B} 1 (x). 



(9) 
(10) 
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These two transformations are referred to [11] as the 
quantum transformation and the background transfor- 
mation, respectively. 

The minimal Landau gauge (in the continuum) is ob- 
tained [12] by considering stationary points of the mini- 
mizing functional 

£[A, g] = J d d x Tr { A™ (x) A™ (x) } . (11) 

Indeed, the first variation with respect to the gauge trans- 
formation g(x) gives 

£[A,g] « £[A, 1] + 2 J d d x Tr{ A^x) D^A] 7 (x) } 

= £[A,l}-2 J d d xTi{ 1 (x)d^(x)} , (12) 

where we used Eq. (5), the relation 

Ty{a,(x) [A^x), 7 (x)]} = (13) 

and integration by parts. (As is usually done, we make 
the assumption that the boundary term in the integration 
by parts gives a null contribution.) Thus, a stationary 
point of the functional (11) satisfies the condition 



Tr{ nd^A^x)} = 



(14) 



which is equivalent to Eq. (2) for A (or) = 0. 

Working in a similar way, one can also obtain the mini- 
mal Landau background gauge. Indeed, the minimization 
of the functional [8] 

£[Q,9] = J d*xTr{Qto{x)QW{x)} , (15) 
yields the variation 

£[Q,g}n£[Q,l] +2 J d d x Tr{ Q^x) D^B] 7 (x) 

+ iQ»(x)[Q fl (x),j(x)]} , (16) 

if we use the gauge transformation (7). The above ex- 
pression may be written as 

E[Q,g] « E[Q,1] - 2 J d d x Tr| 7(2;) D^[B] Q^x) j 

(17) 

if we again integrate by parts, use Eq. (13) and note the 
relation 

Tt{q^x)[B^(x),^)}} = -Tr{ 7 (x)[B^x),Q^x)}} . 

(18) 



Thus, in this case, the stationarity condition implies the 
gauge-fixing relation 



Tr{T b D^[B]Q^(x)} = 



(19) 



which is equivalent to Eq. (3) for A(x) = 0. Clearly, 
for a null background, i.e. B^(x) = 0, the minimizing 
functional (15) coincides with the usual Landau-gauge 
functional (11) and the gauge condition (14) is recovered. 

More in general one should note that, by considering 
quadratic terms in (x) and B^ (x) , there are only three 
terms that can contribute to the minimizing functional of 
the minimal Landau background gauge, i.e. Q M (x) Q^x), 
Qn(x) -B M (x) and B M (x) B^x). However, if one wants to 
obtain the minimal Landau-gauge functional (11) in the 
limit B^ix) — > 0, then the minimizing functional £[Q,g] 
in Eq. (15) is the only choice at our disposal. In this 
sense, the minimizing functional £[Q,g] is unique. More- 
over, of the eight natural sets of gauge transformations 
for the quantum field and the background field (see dis- 
cussion above), one can verify that only the quantum 
transformation (7)-(8) and the set 



Q { "\x) = Q^x) + D^BMx) 



(20) 



B^{x)=B lx {x) +i[Q„(x), j(x)] (21) 

yield the gauge condition (19). Of course, if one lifts 
the requirement of recovering the functional (11) for 
B^{x) = 0, then the minimal background Landau gauge 
can also be implemented by considering for example 
the minimizing functional / d d x Tr{ Q ( ° ] (x) B^ (x) } 
with the gauge transformation Q < jf\x) = Q^{x) and 
B { °\x) = B fl (x) + D^BMx) + 1 [Q At (x), 7 (x)]. 

The above results may be easily extended to the lattice 
formulation of Yang-Mills theories. To this end, we write 
the link variables entering the lattice action as [13] 



We also set 



U»(x) = ^(x)Vv(x) . 

U^(x) = exp[iaA M (x)] 
W M (x) = exp[iaQ^(x)] 
V M (x) = exp[iaB M (ir)] . 



(22) 

(23) 
(24) 
(25) 



where a is the lattice spacing. At the same time, we 
define [14] 



U,(x) t/t(x) 



2ia 



(26) 



traceless 
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and similarly for Q^(x) and B^(x). Then, Eq. (1) is 
immediately recovered, modulo discretization effects. 
The lattice gauge transformation 

Ujf\x) = g(x)U^x)g\x + aeJ (27) 

can also be split among the quantum link W ll (x) and 
the background link V^(x). For example, the quantum 
transformation (7)-(8) is obtained by considering 

Wjf\x) = g(x) W^x) V^x) g\x + ae,) V}(x) (28) 

V^(x)=V,(x), (29) 
I 



in agreement with Eq. (7). 

One can also define a minimizing functional for the 
Landau background gauge on the lattice. Indeed, in the 
limit of small lattice spacing a, the functional 

£[W,g] = -^JiW^i) (34) 

X.fl 

is equivalent to 

£[W,g] « a 2 ]T Tr{Q(f)(x)Qjf)(x) } , (35) 

modulo constant terms. (Here we use 3? to indicate 
the real part.) At the same time, for V^(x) — 1 and 
wjf\x) — Jjjf\x) we recover the usual minimizing func- 
I 



while for the background transformation (9)-(10) we 
have 

Wjf\x)=g(x)W li (x)g\x) (30) 

V^{x)=g{x)V lx {x)g\x + ae ll ). (31) 

Clearly, in both cases the link variable U^{x) transforms 
as in Eq. (27). Moreover, using Eqs. (23)-(25) and the 
lattice definitions of the fields A jt (x). Qf,{x) and B^(x) 
in terms of the link variables U^(x), W^(x) and V^(x), 
one recovers Eqs. (7)-(10) when an infinitesimal gauge 
transformation (6) is considered. For example, Eq. (28) 
gives 

(32) 
(33) 

I 

tional for the Landau-gauge condition [5] 

£{U, g] = ^Tr{g(x)U fl (x)g(x + ae fl )} . (36) 

Also, if wjf^ (x) transforms as in Eq. (28) and we consider 
an infinitesimal gauge transformation (6) we find 

£[W,g]^£[W,l] - i ]T 3Tr{ 7 (x) [U^x)V}(x) 

- v li x - ae n) u ni x - ae v) ] } i (37) 

where 3 indicates the imaginary part. As a consequence, 
a stationary point of the minimizing functional (34) im- 
plies the gauge condition 



Wff\x) « [1 - i 7 (x)} [1 + taQ^x)} [1 + iaB^x)} [1 + i^x + ae^)] [1 - laB^x)] 

« 1 + ia{d^j{x) + Q^x) + i[Q^{x),j{x)] + i [B^x), j{x)] } = 1 + iaQ^\x) , 



= T r{ T b J2 [ W ^ x ) ~ W ^ x ) ~ V K X ~ ae ») M x - ae ») + U K X - ae p) V,(x - ae,) ] } , (38) 

where we used the Hermiticity of the generators T b . Finally, by adding and subtracting Tr{Tb [ W fJ- (x — ae^) — 
W/H 3 - ~~ ae n) } } we find that the null quantity in the above equation can be written conveniently as the sum of two 
terms. The first one is taken as Tr{ T b [W^( x ) ~ W^{x) ~ W^(x— ae^) + W^(x— ae^) ] } and is equal (at leading 
order in the lattice spacing a) to 

2iaTr{T b J2 [QM - Q„(a; - ae„)] } « 2ia 2 Tv{T b ]T 0„Q„(a;)} . (39) 
The second term is then given by 
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Tr{ T b [ U„{x - ae^) Vfa - ae^) + Ufa - ae^) V„(x - ae^) 

-Vfa - ae^U^x - ae^) - V^x - ae M ) Ufa - ae M ) ] | 



(40) 



Note that for a null background field, i.e. B^(x) — 
and Vft(x) — 1, the quantity above is identically 
zero. In this case, we have Q^{x) — A^(x) [i.e. 
W^x) = Uft(x)} and the gauge condition (38) becomes 
[see also Eq. (39)] the usual lattice Landau-gauge con- 
dition Tr{ T h J2„ \M X ) - M x - ae p)l I = °- In thc 
B^{x) ^ case and in the limit of small lattice spacing a, 
one can check that the quantity (40) is, at leading order, 
equal to the expression 



2a 2 Tr{T b V[B M (x), Q^x)]} 



(41) 



Thus, the stationarity condition (38) implies (again at 
leading order in a) 

Tr{T b ^T^Q^x) + i[B^(x),Q^x)}} = 0, (42) 

in agreement with Eq. (19). 

As discussed above, given a fixed lattice configuration 
{Uft(x)}, the usual minimal Landau gauge may be im- 
posed by numerically minimizing the functional (36). In 
particular, by considering local updates for the gauge- 
fixing transformation {g(x)} it is easy to verify that, for 
a given site y, the contribution of g(y) to the minimizing 
functional may be written as [15] 

£[U,g] = constant + 3? Tr j g(y) h(x) j (43) 

with h(x) = E^ix) + Ufa-ae^) + Ufa) + U„{x- 
ae M )]. Then, different gauge-fixing algorithms correspond 
to different choices for the iterative updates of the gauge 
transformation g(y) in Eq. (43). 

In the case of the minimal Landau background gauge, 
one can consider the minimizing functional £ [W, g], de- 
fined in Eqs. (34) and e.g. (28), where {W^(x)} and 
{V^(x)} are given (i.e. fixed) quantum and background 
configurations respectively. It is important to stress that 
also in this case the contribution of g{y) to the minimiz- 
ing functional ^[W 7 , g] may be written as in Eq. (43). In 
this case, the quantity h(x) is equal to 

h ( x ) = [w^) + U t( x ~ ae n) V ^ x ~ ae n) 

+ Wfa) + Vfa - ae^) U^x - ae^) ] . (44) 



I 

Thus, all formulae used for the minimal Landau back- 
ground gauge are natural generalizations of the formulae 
used for the usual minimal Landau gauge. This implies 
that, at least for sufficiently smooth background config- 
urations {V^(x)}, we should expect similar convergence 
of the gauge-fixing algorithms for these two gauge-fixing 
conditions. 

In order to verify this, we have carried out some tests 
in thc SU(2) case, considering lattice volumes V = 8 4 and 
V = 16 4 with a lattice coupling /3 = 2.2, corresponding to 
a lattice spacing a of about 0.210 fermi. This means that 
the thermalized configurations {[/^(x)} are reasonably 
"rough" and provide a good test for the gauge-fixing al- 
gorithm employed. For the background- field {V^(x)} we 
have considered three types of configurations with three 
setups each, namely [here, <jj are the three Pauli matri- 
ces, with o 3 being the diagonal one]: 

a) random center configuration (RCC) V^(x) = 
±1, which can be interpreted as a random configu- 
ration of thin vortices [16], with, on average, 10%, 
30% or 50% of the links equal to -1; 

b) random Abelian configuration (RAC) V fl (x) = 
exp [iQ(x) as], which may be interpreted as a ran- 
dom configuration of Abelian monopoles [17], with 
the angle 9(x) uniformly distributed in the interval 
[0, 27t/] and / equal to 0.1, 0.3 or 0.5; 

c) super-instanton configuration (SIC) [18] given 
by V 2 (x) = exp [i c mm(xi,N — Xi) J2j a j I ^ 3 N 2 ] 
and Vfj,(x) = otherwise, with c = 0.01, 0.05 or 0.1, 
where N is the number of lattice sites per direction. 

For the two lattice volumes above, we consider ten 
gauge-field configurations and, in each case, we fix the 
minimal background Landau gauge, using thc stochastic- 
overrelaxation algorithm [15], for the nine choices of 
background fields described above. The number of mini- 
mizing sweeps necessary to achieve the prescribed accu- 
racy was then compared to that used in the case of a null 
background (i.e. Landau gauge) . Here we stop the gauge- 
fixing algorithm when the average magnitude squared of 
the quantity on the r.h.s. of Eq. (38) is smaller than 
10~ 14 . Note that we tuned the stochastic-overrelaxation 
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B„(a:) 


8 4 


16 4 




aver. 


min. 


max. 


aver. 


min. 


max. 


null background 


217 


190 


290 


508 


396 


773 


RCC 10% 


348 


190 


685 


976 


503 


1729 


RCC 30% 


624 


342 


1391 


1344 


818 


1979 


RCC 50% 


647 


444 


1032 


1711 


1002 


2714 


RAC / = 0.1 


224 


191 


323 


677 


417 


1226 


RAC / = 0.3 


326 


190 


1112 


582 


436 


967 


RAC / = 0.5 


401 


279 


595 


813 


494 


1495 


SIC c = 0.01 


637 


372 


855 


1852 


1238 


3503 


SIC c = 0.05 


188 


172 


256 


520 


344 


808 


SIC c = 0.1 


177 


170 


203 


365 


343 


430 



TABLE I. Average, minumum and maximum number of 
sweeps necessary to achieve the prescribed accuracy for the 
two lattice volumes and for the nine different background 
fields considered in our tests (see description in the text). For 
a comparison, we also include the case of a null background. 

algorithm in the case of a null background, setting the 
parameter p of the algorithm (see [15]) equal to 0.83 for 
V = 8 4 and to 0.91 for V = 16 4 . The same setup was 
then used for non-zero backgrounds. Results of these 
tests are shown in Table I. One sees that the conver- 
gence of the gauge-fixing algorithm for a non-zero back- 



ground is indeed similar to the case of the usual minimal 
Landau gauge. Of course, by tuning the parameter p 
also in the general case, one can improve the results. In 
fact, e.g. for V = 8 4 and background RCC 30%, we find 
that with p — 0.92 the number of sweeps decreases con- 
siderably, being between 418 and 653, with an average 
value of about 460. Similarly, for V = 16 4 and the SIC 
background with c = 0.01, we obtain for p = 0.96 that 
the number of sweeps is between 794 and 1934, with an 
average value of about 1001. 

The above results indicate that numerical simulations 
in the minimal Landau background gauge are indeed fea- 
sible. One should also stress that the extension of the 
method presented here to the case of the minimal covari- 
ant background gauge is, in principle, straightforward 
[6]. This extension, as well as the numerical evaluation 
of Green functions in minimal Landau background gauge, 
is postponed to future studies. 
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